Complex collective states in a one-dimensional two-atom system 
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We consider a pair of identical two-level atoms interacting with a scalar field in one dimension, 
separated by a distance X2i- We restrict our attention to states where one atom is excited and 
the other is in the ground state, in symmetric or anti-symmetric combinations. We obtain exact 
collective decaying states, belonging to a complex spectral representation of the Hamiltonian. The 
imaginary parts of the eigenvalues give the decay rates, and the real parts give the average energy of 
the collective states. In one dimension there is strong interference between the fields emitted by the 
atoms, leading to long-range cooperative effects. The decay rates and the energy oscillate with the 
distance X2i- Depending on X21, the decay rates will either decrease, vanish or increase as compared 
with the one-atom decay rate. We have sub- and super-radiance at periodic intervals. Our model 
may be used to study two-cavity electron wave-guides. The vanishing of the collective decay rates 
then suggests the possibility of obtaining stable configurations, where an electron is trapped inside 
the two cavities. 

PACS numbers: 03.65.-w, 32.80.-t, 73.63.-b 



I. INTRODUCTION 

Systems of interacting atoms form collective states 
where the atoms behave differently from isolated ones 
Q, 0- For atoms in their ground states the collective 
effects are relatively small. They produce the van-der 
Waals or Casimir-Polder forces between atoms Q, asso- 
ciated with the cloud of virtual photons surrounding the 
atoms. 

When the atoms are in excited states they can ex- 
change real photons originated by spontaneous emis- 
sion. Depending on the situation, the real photons 
can give strong collective effects, altering the forces be- 
tween atoms and also their rate of spontaneous emission 
0- El • For example for identical atoms in three dimen- 
sions separated by small distances, the decay rate will 
essentially double or vanish, depending on whether the 
initial state is symmetric or antisymmetric with respect 
to exchange of atoms 0, |(| . We have "super-radiance" or 
"sub-radiance," respectively. If the atoms are not iden- 
tical 0, or if the distance between the atoms is larger 
than their characteristic wavelengths, the effects become 
much smaller 

Most studies on two-atoms (see, e.g., Refs. 0,0, IE El 
IE IE 0> IE IE E3 ) focus on three-dimensional systems. In 
this paper we will consider a one-dimensional system. We 
will show that the exchange of real photons gives strong 
collective effects even for large separations between the 
atoms. 

Our system is analogous to electron wave guid es con- 
sisting of two cavities connected by a lead |l7lll8| . Hence 
the effects we will discuss may be studied experimentally. 

We will consider a simplified model where we have two 
identical two-level atoms, with basis states where one 
atom is excited, while the other is in the ground state. 
We will use the dipole and rotating-wave approximations 
for the interaction with the field. 



We will describe the collective two-atom states through 
complex eigenstates of the Hamiltonian 0, 0, 0, 0] . 
The complex eigenstates decay exponentially in time, 
breaking time-symmetry. The real and imaginary parts 
of the complex eigenvalues give the average energy of the 
collective states and the emission rates, respectively. 

The Hamiltonian, being a Hermitian operator, can 
only have complex eigenvalues if the eigenstates do not 
belong to a Hilbert space. In one-atom systems the non- 
Hilbertian nature of the complex states is manifested in 
their field intensity, which includes a factor growing expo- 
nentially with the distance from the atom. This growth 
in space is related to the exponential decay of the atom 
in time |15| . The exponential field inside the light-cone 
of the atom represents the real emitted photons. As we 
will show, this field has physical effects, which can be 
seen adding a second atom. 

Outside the light-cone, the field associated with the 
complex eigenstates grows exponentially without trunca- 
tion. However, including the complete set of eigenstates 
in the complex representation of the Hamiltonian, the 
exponential field outside the light-cone is cancelled by 
renormalized field states ^E]. This is consistent with 
causality, because the field further away from the atom 
is emitted earlier. Going further away, one reaches the 
point corresponding to the time where the atom was ex- 
cited. At this point the field stops growing. 

One can introduce complex states that are truncated 
outside the light-cone, using distributions dependent on 
the test functions or observables. These are considered 
in Ref. 

In the two-atom system, the real photons emitted by 
each atom are absorbed by the other atom. In d dimen- 
sions the emitted field includes a l/r d_1 decrease factor 
with distance. Hence the field has a strong effect in d = 1 
dimensions, as compared with d — 2 or 3 dimensions. 

We will show that in one dimension the decay rates of 
the collective states oscillate with the distance between 
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the atoms. In contrast to Dicke's states 1] , both symmet- 
ric and antisymmetric states of the two atoms can become 
sub-radiant and super-radiant as the distance between 
the two atoms is varied. For distances that are integer 
multiples of the atom wavelength, the collective decay 
rates vanish, leading to stable collective states. These 
states can trap field energy between the atoms. 

In Sec. [H] we briefly discuss the complex representa- 
tion of the Hamiltonian for a one-atom system. In Sec. 
IIIII we introduce our two-atom model and its complex 
collective eigenstates. In Sees. II VI and IVl we discuss the 
emergence of the collective states and the bouncing of 
photons between the atoms. In Sec. IVII we consider the 
decay rate and average energy of the collective states as 
a function of the distance between the atoms. We dis- 
cuss super-radiance and sub-radiance, including stable 
collective states mentioned above. We also give a heuris- 
tic discussion on the force between the atoms. In Sec. 
IVII I we discuss the mapping of our model to a two-cavity 
electron wave guide, and show that this system allows an 
approximate stable collective state. 



II. ONE-ATOM SYSTEM 

In order to introduce the complex spectral representa- 
tion of the Hamiltonian, we consider first a single two- 
level atom interacting with a field in one-dimensional 
space. This is the Friedrichs-Lee model in one dimen- 
sion. We briefly review the main results. More details 
can be found in Ref . [l4| . 

The Hamiltonian is given by 

H = H +XV (1) 
= Wl |l)<l| + ^^|fc)(fe|+A^V4(|fc)(l| + 



L — > oo. For L finite the momenta k are discrete. In the 
limit L — > oo they become continuous, i.e., 



2tt 



E 



dk. 



(2) 



We have (a\b) = 5 a b = Kroenecker delta. In the limit 

L — > oo, 



^-5 k ,w - S(k - k'). 
Ztt 



(3) 



The interaction term is obtained through the dipole 
approximation as well as the rotating-wave approxima- 
tion. The potential V4 is of order L^ 1 / 2 . For convenience 
we write 



V k = {2n/L) 1 / 2 v k , 



(4) 



where v k is of order 1 in the continuous spectrum limit 
L — > oo. As a specific example we will assume that [l9L 

M 



Vk 



v(uj k ) = 



1/2 



[1 + (wt/wjf) 2 ]" 



(5) 



with n = 1. The constant ujJj determines the range of 
the interaction. We shall assume that the interaction is 
of short range, i.e., u>m ^ <*>\. 
The state |1) is unstable if 



dk^ 



(6) 



Otherwise, it is stable [23. Hereafter we will consider 
the unstable case. 

In the unstable case one can construct renormalized 
field eigenstates that diagonalize the Hamiltonian as 



where we put c = % = 1. The state |1) represents the 
bare atom in its excited level with no field present, while 
the state \k) represents a bare field mode ("photon") of 
momentum k together with the atom in its ground state 
(see Figure P). 



A/VW 



H 



where 



k 



hm|^) = |fe) 



(7) 



(8) 



and hereafter we use the summation over field modes in 
the sense of Eq. @ . The index ± refers to either "in" or 
"out" scattering eigenstates. 
The explicit form of the eigenstates is given by [l4| 



II > lk> 



FIG. 1: One-atom system 



The energy of the ground state is chosen to be zero; 
u)\ is the bare energy of the excited level and u>k = |fc| 
is the photon energy. The coupling constant A <C 1 is 
dimensionless. We assume periodic boundary conditions. 
We put the system in a "box" of size L and take the limit 



l^> = I*) + 



E 



1 



u>k — u>i ± it 



-10 



(9) 



where e is an infinitesimal positive number. The limit 
e — * is taken after the limit L — > 00. In Eq. ©, 



77 ± (cj) = uj — u>i — 



A 2 ^ 



■f (u> - UJk'Y 



uj — <jj\ 



2 1 dk' 

/o ( w 



A V k> 



(10) 
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is the inverse of Green's function. The + (or — ) super- 
script in Eq. I|1U|) indicates analytic continuation from the 
upper (or lower) half- plane of u 14]. Using the complex 
delta-function 5c [Tl| we can write 



f (w - k' 



zblmcj > 



(oj - fe')± I ( w - T 2iriS c (k' - w), ±Imw < 



(11) 



When a; is real, we have 



dk'- 



(12) 



With our form factor and small A, Green's function has 
one pole Z\ in the lower half plane, i.e., r/ + (zi) = for 



(13) 



The negative imaginary part —171 describes decay for t > 
0. The real part u)\ gives the shifted average energy of the 
excited state. For the other branch we have 77" (z*) = 0, 
with z\ describing decay for t < 0. 

Note that in the representation {JJ the decay rate and 
shifted energy of the excited state do not appear in the 
spectrum. One can incorporate z\ (or z*) into the spec- 
trum by extracting the residue of Eq. (0 at the pole 
Wfc = z\ (or zf). This gives the complex spectral decom- 
positions in til mni: 



H 



T l/2 



ID 



k 
k 

1 1 an 

m 



(14) 



The state \4>i) and its dual (01 1 are complex eigenstates 
of H . Their explicit forms are 



|^i) = [iV 1 *] 1 / 2 [|l)+^|fc) 



(Zl - LUkY 

AT4 



(15) 



(z* -u k y 



The state \4>k) h as the same form as the state 
the replacement 



1 



LO k - Zl 



V + (u k ) 



with 



(16) 



and similarly, the state \4> k ) has the same form as the 
state \4>k) with the complex conjugate replacement. 

The states in Eq. I|14f) form a bi-orthormal set, with 
the relations 



1, 







#14 



> = 5k,k> 



(17) 



and their complex-conjugate relations. 



III. TWO-ATOM SYSTEM 

In this Section we discuss the complex spectral repre- 
sentation of a two-atom system with Hamiltonian 

H = wi|l)(l|+W2|2)(2|+5^« fc |fc)<*| 

fe 

+ ^2XiV k (e ikx ^\l)(k\+e- ik ^\k)(l\) 

k 

+ Y. x z v k( elkX2 \ 2 )( k \ + e ~ lkX2 \ k )( 2 \) ( 18 ) 

The state |1) represents atom 1 in its excited state, while 
atom 2 is in the ground state and no held is present. 
Conversely, the state |2) represents atom 2 in its excited 
state, while atom 1 is in the ground state and no field 
is present. The state represents a field mode k with 
both atom 1 and atom 2 are in their ground states (see 
Figure|2J). The atoms 1 and 2 are located at the positions 
X\ and x<ii respectively. We use the potential in Eq. JSJ. 



Il> 



12 > 



AA/W 



lk> 



FIG. 2: Two-atom system 

We will first assume that the two atoms are at fixed 
positions, so that the distance between them 



X21 = \x 2 - xi\ 



(19) 



is fixed. This can happen if the atoms are heavy. A 
system of two fixed atoms is analogous to a two-cavity 
waveguide (see Sec. IVIIl! . 

We will consider the case where the two atoms are 
identical, 



Ai — A2 — A 



(20) 



We introduce the symmetric and antisymmetric states 

(21) 



n |1) + |2) I v |1)-|2) 



y/2 ' 



V2 
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which are eigenstates of the unperturbed Hamiltonian 
H = H(X = 0) as 



Ho\j) = j = s,a 

UJg = UJ a = UJl. 

We will use the notation 

_ J 1, for j = s 
°3 -\ -1, for j = a, 



(22) 
(23) 



(24) 



With this notation we have \j) = (|1) + e>j |2))/\/2. 

As for the one-atom system, we can diagonalize the 
Hamiltonian as 



(25) 



where 



\Fk) = \k)+(^\s)+^ k \a)+Y,^ k \k') (26) 



and 



8± = — XVk ( c ikxi I a c ikx A C>7) 



(28) 



»7?M 



= w - u>x - 2 / dk'- (1 + cjj cos(fc'a;2i)) 

(29) 

Following a procedure similar to one found in Ref. [l4| . 
one can show that the new diagonalized states \Fj^) sat- 
isfy the orthogonality and completeness relations 

E I^X^I = |1><1| + I2><2| + £ (30) 

fe k 

{*t\F£) = {K\Fv)=*KV ( 31 ) 

Green's function [?7^~(w)] _1 has poles in the lower half- 
plane, and conversely, [f^~(w)] _1 has poles in the upper 
half-plane. From now on we discuss only the + branch 
with poles on the lower half-plane. 

A new feature with respect to the one-atom system is 
that due to the cosine term in Eq. (|29() . there are many 
poles of Green's function, as shown in Fig. [31 We label 
the poles as 



(32) 



where n is an integer. 



The poles Zj, n are solutions of the equation 

Vl(zj,n) = (33) 

In the following we discuss this equation and its solutions. 
From Eqs. lJTTJl and we obtain 



UJl + 2 



dk 



X 2 i, 2 



- (1 + a-j cos(k'x 2 i)) 



4Tri\ 2 [v Zj n ] 2 (1 + (jj cos{zj^ n x 2 i)) 



(34) 



Note that cosine in the last term includes the factor 
exp(7j.„cc2i), which grows exponentially with the dis- 
tance X21 between the atoms. Assuming weak coupling 
and taking only the pole contribution in the k! integral 
we obtain the set of equations 

G)j >n + 2Tr[\v(ujj, n )] 2 <jje 1j ' nX21 sin(uij >n X2i) (35) 



7j,„ w 2ir[\v(LU : j, n )} 2 [1 + <7,,-e 7i ' nX21 cos(Qj, n x 2 i)] (36) 

In Fig. |21 the pole of [^(uj)]^ 1 with real part closest to 
the unperturbed frequency u>\ is also closest to the real 
axis. We call this pole z s ^ = z s . Similarly, for [77+ (w)] _1 
we denote the pole closest to u>i as z a> o — z a . Both these 
poles are obtained by a perturbation expansion around 
A = 0. We have 



2.7,0 



Zj — > u)\ as A — > 0. 



(37) 



If X21 is not too large {x%\ ~ r y 1 1 or smaller), one can 
show that the poles Zj n are given by 



Zj + 2mr/x 2 i + 5zj <n , 

zj + (2n + <Tj)n/x 2 i + 8zj <n , 



for (jjn > 
for UjU < 0. 



(38) 



where S Zj'n IS £111 0(A 2 ) correction. The approximate 
value Re(Zj n — z j) predicted by this equation agrees with 
Fig. |3 (for j = s) and a similar figure for j = a, which 
we omit. 

We write the poles Zj as 



(39) 



The poles Zj , having the smallest decay rates jj , will give 
a dominant contribution to the time evolution after a few 
bounces of the field between the atoms. In this way, the 
complex collective states defined in Eq. l|40(l emerge. 

As in the one-atom system we can obtain complex 
eigenstates \<f)j) of the total Hamiltonian, such that 



H\^)=z j \^j) 
Their explicit forms are given by 

\4>i) 



(40) 



N 



1/2 



li) 



2~ 1/2 XV k , _ tkxi 



-ikx2 



\k) 



(41) 



5 




12 3 4 



FIG. 3: Contour plot of log(l/\r/f (z)\). The x and y axes 
are Re(z) and Im(z), respectively. The contours concentrate 
around the poles z s ,n of Green's function. Parameters are 
X21 = 29.025, wi = 2, A = 0.05, ui M = 5. 



1 2 3 4 5 6 7 



FIG. 4: Log plot of the survival probability Pi(t) (solid line) 
and complex collective state component Pi, zs (t) (dashed line). 
The crosses indicate the decay rate 2y sl , *y sl — 0.0233, be- 
tween t = and t — 1 (see Sec. 0. Time is in units of 
221- The distance between atoms is X21 = 29.025. Other 
parameters are u>i — 2, A = 0.05, ujm = 5, and L = 500. 



where 



N: 



f [(Zj - uj k )+} 



;(1 + <Tj cosfcz^i) 



(42) 



For these states we have \<f>j) — > \j) as A — > 0. 

The dual states satisfying (<j>j\H = Zj(<j>j\ are given by 



1/2 



01 +E 



2- 1/2 AT4 , ikxi 



cr.,e 



tkx-2 



(k\ 

(43) 



Like in the one-atom case, we have the complex spectral 
representation 



Ei^ + >^ 



(44) 



where \F k ) has the same form as the state \Fj[) with the 



replacement 



LUk — Zj 



(45) 



We have as well the complex-conjugate representation, 
taking the complex-conjugates of Eqs. (|4^l) . l|43|) . 



IV. EMERGENCE OF THE COMPLEX 
COLLECTIVE STATES 

Time evolution of the two atom system can be solved 
by using Eq. I|25|l or Eq. I|44|) . As an example we assume 
the atoms are initially in the symmetric state \s) and 



the initial field is zero (similar calculations can be done 
if the initial state is \a}). We will calculate the survival 
probability of state |1), 



Pi(t) 



le 



lHt \s)\ 2 



(46) 



Before we go into details, we can guess the behavior the 
system will show. Since the initial state is symmetric, 
the following discussion also applies with atoms 1 and 2 
exchanged. Say atom 1 is to the left of atom 2. At the 
beginning, atom 1 decays and emits a field. Half of this 
field will be radiated away to the left, while the other half 
will reach and excite atom 2. Atom 2 will then decay and 
emit its own field, part of which will be radiated away to 
the right, the rest going to the left, back towards atom 1. 
Continuing this process, we see that energy will bounce 
back and forth between the two atoms. As time passes, 
this energy will decrease due to the outgoing radiation. 
Eventually both atoms will decay to the ground state. 
Noting that the time it takes for the field of one atom 
to reach the other atom is t — X21 (with c = 1) we 
conclude that, as it decreases, the survival probability 
should oscillate with period X2i- 

This behavior is shown in Fig. This was obtained 
through a numerical solution of Schrodinger's equation. 
The field was discretized into 2501 modes. The eigenval- 
ues and eigenfunctions of the Hamiltonian matrix were 
obtained using tri-diagonalization and the "QL" method 
[23| . This allowed us to calculate explicitly the operator 
exp(-iHt). For this and the subsequent numerical plots 
we used the following parameters: oji = 2, A = 0.05, 
wjk = 5. Other parameters are indicated in each figure. 

In order to calculate the survival probability we start 
with Eq. (|23|) to obtain 



Pi(t) 



(He 



-iHt 



Ei^ + n^ + i s > 



6 



(i|^ + X^ + k) 



(47) 



where used the fact that odd functions of k vanish under 
the summation. For later use we define the amplitude in 
Eq. |g7J| as 



~™ kt +, V \,, (l + cosfcs 2l ) (48) 



The dominant contribution to P\ (t) will come from the 
poles of Green's function, shown in Fig. [3J The different 
pole contributions should add up to give the bounces 
seen in Fig. But rather than computing all the pole 
contributions, we will follow an easier method in Sec. Ivl 

Here we will focus on the pole z s . As mentioned be- 
fore, this will give the dominant contribution after some 
bounces, since it gives the slowest decay rate. It is this 
pole contribution that is extracted in the representa- 
tion H44[l. Using this representation, and noting that 
(<Pa\s) — we have 



energy and the field P(x,t) around the atoms starts to 
approach the field intensity due to the collective state 
given by 



P zs (x,t) = \{ip(x)\(j) s ) exp(-iz s t)(<f>s\s) 



(53) 



(see Fig. 0). The collective state decays exponentially 
(Fig. 01. 



'Ms 



-iz s t i 



J2(l\F+)e-^{F+\t 



(49) 



The second term contains contributions from the poles 
other than z s of Green's function [77+ (w)] -1 as well as 
contributions coming from the branch cut of this func- 
tion. Neglecting all these contributions we obtain 



Pi(t) & P 1<za (t) = (l|0 a )e 



(50) 



This is represented by the dashed line in Fig. After 
a few bounces the initial state \s) reaches the collective 
state \(f> s ). 

We turn to the time evolution of the field. Defining 
the state 



1 



- ikx 



(2w fc L)V2 



«|fc), 



(51) 



the intensity of the field in space-time can be written as 

P(x,t) = mx)\e- im \s)\ 2 (52) 

Again we calculated this using the numerical solution 
of Schrodinger's equation. The intensity of the field is 
plotted in Figs. I5l7l for different times. At the beginning, 
both atoms emit their fields spontaneously. Each field 
has an exponentially growing envelope (plus corrections 
due to the initial dressing processes (l5|). which stops at 
the light cone \x — Xi\ = t (Fig. [5). 

After each emitted field reaches the neighbor atom, ab- 
sorption and re-emission occur. The two atoms exchange 



FIG. 5: Field intensity P(x,t) for t = 0.32 X21. The atoms 
are located at x\ = and X2 = 1. Space coordinate x is in 
units of £21 and P(x,t) is dimensionless. The parameters are 
the same as in Fig. [1] 




FIG. 6: Field intensity P[x,t) for t = 1.12 X21 (solid line) 
and the complex collective state component P za (x,t) (dashed 
line). Space coordinate x is in units of X21 and P(x,t) is 
dimensionless. Parameters are the same as in Fig. 0] 



In summary, the atoms emit a field growing exponen- 
tially with the distance from them, within their light- 
cones. After the field emitted from each atom reaches 
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FIG. 7: Field intensity P(x,t) for t = 4.02 x 2 i (solid line) 
and the complex collective state component P zs (x,t) (dashed 
line). The outer smaller peaks of P(x, t) come from the initial 
one-atom emission. The inner, larger peaks come from the 
emission after the first exchange of energy between the atoms. 
P{x,i) asymptotically approaches P zs {x,t) as x approaches 
the atoms. They coincide in the region between the atoms 
(between x = and x = 1). Parameters are the same as in 
Fig. H 



the other one, the collective state with complex energy 
Zj emerges. 

As we discuss now, the exponential field has a strong 
influence on zj. The field amplitude associated with 
the collective states is given by (k\<fcj). This amplitude 
in turn determines Zj through its interaction with the 
atoms. We have 




ikx± 



ikx2~\ 



(/#,) (54) 



where we used Eqs. <j4"5|) . (|4*TJ) . and (J2HJ, with lu = zj. 
Since (k\(f>j) are functions of Zj, this is a self-consistent 
relation. The exponential component of the field is seen 
in the approximate equations (|33|) and for zj t0 — Zj , 
which include the factor exp (7^x21). Due to the expo- 
nential nature of this factor, the pole Zj may deviate 
substantially from the one-atom pole Z\, 

In spite of the exponential factor, for increasing X21 
the equations lj3*BTl and l|3*B|) can still have solutions if jj 
decreases as 



(55) 



for large £21 • The decrease of jj with increasing X21 is 
seen in Fig. |S] 



V. BOUNCES 

In this section we describe the energy bounces between 
the atoms, seen in the survival probability of each atom. 



As shown in Fig. the decay rate of Pi(t) changes 
abruptly at t — x%\. For t > X21 the decay rate quickly 
approaches the collective decay rate j s - The wiggling 
of the decay rate shows the absorption and re-emission 
of the fields, or in other words the energy bounces. For 
< t < X21 the decay rate should be close to the one- 
particle decay rate. 

To analyze the energy bounces and the decay for < 
t < X21 we first note that 

vfi k ) - Vs(k) = 4:TTi\ 2 vl(l + cosfcE 2 i), (56) 

Hence in Eq. Q48[l we can write 

r°° X 2 v 2 
I(t) = 2 / dke-™** k (l + cosfcxzi) 

Jo Wwl 



1 

2tt7 



dk 



1 



Vs(k) vt(k) 



Since e~ tkt vanishes in the lower infinite semi-circle of 
complex k plane for t > 0, we can take the pole contri- 
butions extending the k integration from — 00 to 00 and 
closing the contour with this semi-circle. Only [^(k)}^ 1 
has poles in the lower half-plane. We write rj^(k) as 



r]+(k) = k-u) X - 2 / dk 
Jo 

= v +(k)-A(k) 
where ^(fc) is defined by 
Vf^k) 



,\ 2 vl,(l + COsfc'x 2 l) 

k — k' + ie 



(58) 



k - ui - 2 / dk' 

J k — k' + ie 

„,\ 2 vle lk ' X21 



rr-a + k 21 ) 



and 



A(fc) = ~2iri\ 2 vle lkX21 



(59) 



(60) 



Unlike [r]^(k)] 1 , [f?^(fc)] 1 has only one pole in the 
lower half plane. Let this pole be 



Z s i = UJ s i - l^fsi 



(61) 



This is essentially the pole of one-atom Green's function, 
modified by the overlap of the atomic clouds at the dis- 
tance X2i- For X21 3> Wj" 1 we have 



Z S 1 ~ Z\. 



(62) 



We expand l/r]f(k) as 
1 



V t(k) !,+ (*:) -A(fc) 
A(k) n 



E 



o (V7 1 (k)r +1 



(63) 
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The expansion is possible since 
+27rA 2 ^ sinfcx 2 i + 2niX 2 vl 



o . 1 



2 / dk'V 
Jo 

2„,2 „• ; — , n_- \2„,2 



A 2 t^,(l + COS k'Xil) 



k~k> 



\A(k)\=\lm( V +(k))\<\( V t 1 (k))\. 
Using Eq. Ip5j). Eq. JSJJ is written as 

„— ikt 



(64) 
(65) 



i r 



dk- 



2m 



£ 

n=0 



% (fc) 

(-27ra 2 ^)™e- lfc (*- n:E21 ) 



n+l 



(66) 



In Eq. (|66|) . the pole contributions come from 
l/( 7 ?si(^))™ +1 - For n = 0, e.~ lkt / ^{k) has a simple pole 
in the lower half plane at k = Its effect appears for 
t > 0, when we can close the integration contour in the 
lower half plane. For n = 1, e~ lk ^~ X21 ^ / (^(k)) 2 has a 
double pole. Its effect appears for t > x 2 i- In general, for 
each X21 time step there appears a new pole effect which 
is smaller by A 2 order than the previous pole effect. In 
this way we can explain the wiggling decay rate (Figure 
0. 

As we discuss now, this description of the bounces is 
connected to emergence of the collective state. Approxi- 
mating (for A <C 1) 



Vfiik) ~k- z sl 
the pole contributions in Eq. (|66|) are given by 

1 



(67) 



Io(t) 



-1 
27ri 



1 OO 

— y 

2ni ^ 



dk 



k — z s \ — A(fc) 



dk 



A(ky 



ikt 



-ikt 



(68) 



Taking the residues at the pole k 
expression of the form 



(fc-Z sl )«+! 

z s i we obtain an 



- nxix)f n (t) 



where 



/«(*) 



1 d n 



i\ dk n 



A{k) n e 



n^—ikt 



(69) 



(70) 



We have /„ ~ A 2 ™. Note that the sum stops at n such 
that t < 71X21- Since for weak coupling the terms f n {t) 
become smaller as n increases, after a few bounces we 
have 



where 



Io(t) « I (t) 

OO 



(71) 



(72) 




FIG. 8: The decay rates 7 3 (x) and 7 a (+) oscillating as a 
function of X21 ■ The solid line is the one-atom decay rate 71 = 
0.0235. 7s vanishes for distances close to (2n — 1)tt/uji, and 
7 a for distances close to 2mr/Cji with n integer. An example 
is £21 = 12.7 ~ 8tv/Qi, where *y a vanishes, while y a is large. 
For large X21 the decay rates decrease. 



As shown in Appendix^] we have 
J (t)=iV s e-^ t 
for all i > 0, where 

N, 



l-dA{k)/dk 



k—Z s 



(73) 



(74) 



for weak coupling. Eq. I|73|) shows that the sum of all 
bounces gives the contribution from the collective state 
\<p s ) with eigenvalue z s . 

Eq. (|73|l is consistent with z s giving the slowest expo- 
nential decay. To see this we use Eq. Ij58(l to write the 
equation for z s ^ n as 



or 



1ltl( z s,n) ~ A(z s ,„) = 



Z s l + A(z s ,n) 



(75) 



(76) 



[for X21 > 1 we have z s i ~ z\ and we recover Eqs. 
and for n = a}. 
The function k — z sl — A(fe) in Eq. (|68|) has zeroes 
at k — z s>n . For t — » 00 only the residue at the pole 
k = z s ,o = z s remains and we get 

lim A7V^%(i) = lim Ar s -V z »*Jo(i) = 1 (77) 

t — >oo t — >oo 

which is consistent with Eq. 173|) . 



VI. DECAY RATE AND ENERGY VS. 
DISTANCE 

In this section we investigate the behavior of the com- 
plex eigenvalues of the Hamiltonian Zj for different values 

of X 2 1- 
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FIG. 9: The energies ui s (x) and u) a (+) as a function of £21- 
They oscillate with X2i- The solid line is the one-atom energy 
Qi = 1.985. 



FIG. 10: Survival probability Pi(t) of atom 1 for symmet- 
ric (dashed) and antisymmetric (solid) initial conditions and 
X21 = 12.7. The symmetric state gives rise to a super-radiant, 
stationary collective state. The antisymmetric state gives rise 
to a sub-radiant (stationary) collective state. Before t — 221, 
Pi (t) has the one-atom decay rate. Time t is in units of £21 • 
Pi(t) is dimensionless. 



The equation rQ{z) — can be solved numerically by 
iterations of z — z — iy(z). The imaginary and real parts 
of z thus obtained are shown in Figs. OH and 0D (we used 
the same parameters as in the previous figures). The 
numerical iteration was started around z = lu± so, with 
the exception of two isolated points seen in Fig. 00 the 
solutions obtained are the collective eigenvalues z = Zj. 
Gaps in the graphs are points missed by the numerical 
solution. 

As we see, 7^ and ujj oscillate with X2i- The oscillation 
period is approximately 2it/uji where Cj\ is the one-atom 
renormalized frequency (see Appendix [BJ). 

Due to the oscillations, the collective decay rate can 



FIG. 11: Field intensity P(x,t) for the atoms in a stationary 
collective state. The field between the two atoms, located at 
x — and x = 1 remains trapped, x is in units of X21 ~ 12.7 
and P(x, t) is dimensionless. The initial condition is \a). Time 
is t — 7.09 £21- The wave packets on each side represent the 
field emitted before the atoms formed the collective state. 
After this, there is no emission (we have sub-radiance). 
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0.007 
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FIG. 12: Field intensity P(x,t) for initial condition |s) at 
time t = 7.09 aai- x is in units of 2:21 = 12.7 and P(x,t) is 
dimensionless. The collective state has decayed. The smaller 
peaks on the far sides are the field emitted individually by 
each atom before attaining the collective state. The larger 
peaks correspond to the two-atom collective emission (super- 
radiance) . 



become smaller or larger than the one-atom decay rate 
(solid line in Fig. |SJ). We have sub-radiance and super- 
radiance, respectively. In particular, it is noticeable that 
there are distances at which the decay rates jj vanish 
(see Appendix [Bj . This means that for these distances 
there is no outgoing radiation. The outgoing emitted 
fields of the atoms cancel by destructive interference and 
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FIG. 13: The graph of T a = —dCj s /dx2i (black line) and the 
decay rate — 7 S (dotted line). We see that when — 7 S is a local 
minimum T s has also a local minimum value. 



attractive force between two atoms becomes locally max- 
imum when the collective decay rate is locally maximum. 
The atoms tend to attract each other when they emit the 
field outwards, and tend to repel when the field remains 
trapped between them. 

Also we see in Fig. ^| that there are points x^x for 
which T a vanishes. If dT s j dx2\ < at these points, then 
any small displacement AX21 around x^i creates a force 
in the opposite direction. Thus in this CcLSG RFC stable 
points [if dT s fdxi\ > the points are unstable]. The ex- 
istence of stable points suggests the possibility of having 
a one-dimensional "molecule." This molecule would have 
a lifetime of the order of 7" 1 . 

VII. TWO-CAVITY WAVE GUIDES 



a standing field is trapped between two atoms, storing 
energy. Note that both symmetric and antisymmetric 
initial conditions can give rise to either sub-radiant or 
super-radiant states, since both 7 S and j a oscillate with 

X21- 

The oscillations of the decay rate and the energy shown 
in Figs. [S] and [5] are a unique feature of one-dimensional 
systems. For two or three dimensions, these quantities 
can change significantly only for short distances between 
atoms (see Appendix 0. 

As an example of sub-radiance and super-radiance we 
show numerical simulations with the same parameters 
used in the previous examples, except we choose L = 250 
(to have higher space resolution) and X21 = 12.7. For this 
value of X21, the decay rate of the antisymmetric state 
vanishes while the decay rate of the symmetric state is 
maximum, (see Fig. [SJl. In Fig. ^|we show the survival 
probability of atom 1 for the antisymmetric and sym- 
metric initial conditions, showing the appearance of sta- 
tionary sub-radiant collective state and a super-radiant 
collective state. In Figs. ^2 an d El we show the corre- 
sponding fields. 

We turn to the force between the atoms. Here we will 
only give a heuristic discussion. A more detailed analysis 
requires including the Casimir-Polder or van der Waals 
forces between the atoms, as well as the inertia of the 
atoms, which we are not considering in this paper. 

Since the atoms are unstable, the force between them 
should be time-dependent [24[ . We expect the force will 
decay exponentially during the time scales where the 
collective-state components dominate. For the depen- 
dence on X21 of the force, the quantity 

J-j = —dujj/dx2i (78) 

can give an indication because ujj is the average energy 
of the collective state. 

As we can see in Figure ED oscillates with X21 {J- a 
has a similar behavior). T s > corresponds to a re- 
pulsive force, and T s < to an attractive force. The 
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FIG. 14: A two-cavity waveguide 



As shown in Rcfs. 0,0], a Hamiltonian of the form 
()18J) can be used to describe a two-dimensional electron 
wave guide as seen in Fig. 1141 

This wave guide can be constructed by superposing 
two closed identical cavities and a lead. This forms the 
unperturbed system. The interaction appears as the cav- 
ities and the lead are connected. 

In suitable units the horizontal dimension of the cavi- 
ties is 1, and the vertical dimension is D. The lead has 
a horizontal dimension L — > 00 and a vertical dimension 
W. We consider a non-relativistic electron, neglecting 
the spin. 

If the electron is placed inside a closed cavity, its wave 
functions correspond to discrete cavity modes. The cav- 
ity modes can be labeled as |m,n), where m, n are posi- 
tive integers representing the horizontal and vertical wave 
numbers. The corresponding energies are 

C" = m 2 +n 2 /D 2 (79) 

An electron placed inside the lead (with no cavities) 
has modes that can be labeled as \k, 1} where kL/n is the 
horizontal wave number and I the vertical wave number. 
The energies are 

E k .i = k 2 /ir 2 + l 2 /W 2 (80) 

As L — > 00, k becomes a continuous variable. On the 
other hand, m, n and I are always integers. 

We consider an electron with low energy narrowly cen- 
tered around 

£0 _ £m ,no j-g-^ 
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We assume that 



Eo,i < C < E k ,i 



(82) 



for I > 1 and all k. The electron may propagate through 
the first mode of the lead, but not through the I > 1 
modes. 

We also assume that there are no other cavity modes 
with energy between 25 0)1 and Co- Under these conditions, 
the cavity mode £° behaves essentially like the excited 
state in the Friedrichs-Lee model. It will decay with a 
finite lifetime. This means than an electron inside the 
cavity will escape through the lead. 

The following approximate Hamiltonian is obtained 

m 

Hwg = e°[|l)<l| + |2)<2|]+5^ J B fc ,,|Vfc,i)<V'fc,i| 

k,l 

+ J2Ki (|l>^Me ifcxi + \2)(<<p k , l \e ik *>)+h.c. 

k,l 

(83) 

Here, h.c. means Hermitian conjugate. The cavities are 
centered at x = x\ and x = x%, where x is the horizontal 
coordinate. The states \i) represent the electron inside 
cavity i = 1 or 2, occupying the mode Co- The states 
\ipk,i) are modified lead modes; they essentially represent 
the electron inside the part of the lead that does not 
overlap with the cavities. The terms V fc ° ; represent the 
amplitude of a transition of the electron from this part 
of the lead to the cavities or vice versa. Their detailed 
expression is given in Ref. |l8j . 

Except for the additional index I and the dispersion re- 
lation H8l)f) . which is different from u>k = \k\, the Hamilto- 
nian Hwg is the same as our two-atom Hamiltonian (11811 . 
As the two cavities are identical, we have a system anal- 
ogous to the two identical atoms. Since there is only one 
continuous variable k describing the propagation along 
the lead, we can think of the wave guide system as a 
one-dimensional system, with internal degree of freedom 
I (note that / is discrete). 

Using the results of the Sec. IHII we obtain the equation 
for the complex energy of the collective state 



c° 



oo 

\ " 



-(1 + <7j COSfcX2l) 



(84) 



where = (L/irjV^. We will show that there is a 
solution with vanishing decay rate, corresponding to a 
stable collective state. We follow the procedure shown 
in Appendix [5] For a vanishing decay rate we write 
= I-® — ie, where e > is infinitesimal. This gives 
the following condition on X21 '■ 



1 + <7j COS k()X21 = 

where ko is a wave vector that satisfies 



Ek ,i = Cj, I = 1 



(85) 



(86) 



Through these two equations, X21 becomes a function of 
C° 



X21 = 



(87) 



where n = odd integer for j = s, and n — even integer 
for j = a. The renormalized energy is then given by 
the solution of the integral equation 



I? = Co 



1=1 



to 



E, 



1 + <Ti (x»[hg(B) 



where V means principal part. Similar to Eq. (|B11I) . this 
equation has a solution if the condition 



/•oc x 

e-E 0A >2 

J ° 1=1 



Ek,i — Eqi 



(89) 



is satisfied. If the two cavities are not too close, we can 
replace the interaction v® l by the interaction in a single- 
cavity system. Then, Eq. I|89|l is essentially the condition 
that the electron in an individual cavity has enough en- 
ergy £° to escape through the lead. This condition is 
analogous to Eq. @. 

In summary, adjusting the distance between the cavi- 
ties, so that Eq. ||%7I) is satisfied, we obtain a collective 
stable state where the electron remains trapped inside 
the two cavities, in either a symmetric or antisymmet- 
ric state. The electron is trapped even though it would 
escape if there was only one cavity. 

To obtain this result we neglected the influence of cav- 
ity modes other than £°. The existence of stable configu- 
rations in the wave guide could be verified by other meth- 
ods, including numerical simulations or experiments. 



VIII. CONCLUDING REMARKS 

We have analyzed a two-atom system using complex 
collective eigenstates of the Hamiltonian. Our main 
result is the description of long-range effects in one- 
dimensional space, such as the vanishing of the decay 
rate at regular intervals of distance £21 ■ Another result 
is the application of this model to two-dimensional elec- 
tron wave guides, where a two-cavity configuration can 
be tuned to act as an electron trap. 

In our two-atom model we neglected virtual transi- 
tions. This corresponds to a rotating wave approxima- 
tion. Phenomena such as the existence of collective stable 
states in one-dimensional atoms deserve further study, 
with the inclusion of virtual transitions. Electron wave 
guides, on the other hand, are already well described by 
the type of Hamiltonian we considered, without virtual 
transitions. The description improves if we include more 
cavity modes 0] . 
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Emitted photons are described by an exponentially 
growing field, truncated at the light cone of the atoms. 
This field plays an important role in the two-atom sys- 
tem, giving a strong influence on the lifetime or average 
energy of the collective states. This field is directly re- 
lated to the exponential decay of unstable states, which 
can be regarded as one of the simplest dissipative phe- 
nomena on a microscopic scale. So, in a sense, the forma- 
tion of collective states is a microscopic non-equilibrium 
process, driven by dissipation. 
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APPENDIX A: PROOF OF EQ. fTSl) 

We start with the expression (see Eq. J72J) 

A(k) n 



where 



r 

'<>(<) = / dk 



-ikt 



(Al) 



where C is a clockwise contour surrounding k = z sl . 
Taking the residues at this point we get 



k—z s i 



(A2) 



This is a perturbation expansion around z s \, so it will 
correspond only the contribution from the pole z s and 
not the poles z SlTl . 
We will show that 



d - 

— Jo(t) = -iz s I (t). 



(A3) 



Together with Eq. (J77J) this will prove Eq. (|73J). Starting 
with Eq. (|A2|) and using 



dk- 
we obtain 



-ab = y 



" Qn—l 



1=0 



ll(n-l)\ 



Q k n-l 



A 



Ql 

dk 1 



D 



—I (t) = -iz s Jo(t) - i~h{t) 



(A4) 



(A5) 



-1 d 



Lm = E ^d^ [A{k)n+me ~ M]k ^ (A6) 



Using Eq. (|A4|) again we get 



Lit) = E ^m^iAikr}^ (A?) 



With Eq. (|76[) for z Si q — z s and the Taylor expansion of 
A(z s ) 1 around z s \ we find that the solution of this system 
of equations is 

/,(*) = A(z s ) l I (t) (A8) 
which combined with Eq. I|A5J| proves Eq. (j A3|l . 



APPENDIX B: OSCILLATIONS OF 7j AND uj 3 
WITH X21 

In this Appendix we will show that the decay rates jj 
and energies Cjj of the collective states \4>j) oscillate with 
the distance x 2 i between the atoms. First we will show 
that j s vanishes (comes infinitesimally close to zero) for 
distances 



[^2l] 7s =0 = 

where n is an integer, and 



(2n+ 1)tt 



(Bl) 



u>° = p s ] 7s =o (B2) 
Similarly we will show that decay rate j a vanishes for 



2mr 

F2i] 7 a=o = 77^- 



where 



= Pa] 7a =0- 



(B3) 



(B4) 



We start with the equation 77+ (z s ) = or 

Z s = LUl + 2 



\ 2 v 2 

dk- 1— (l + cosfca; 2 i) (B5) 

[ z s — k) 



Assuming z s — uj° — ie with infinitesimal e we have 
u>° = wi + 2 



\ 2 v 2 

dk— — — (1 + cos kx2i) 



= oji + 2 



uj° — k + ie 
V 



dk\ 2 v 2 k 



X (l + COS kX2l) 

where we used the relation 
l V 



U! + It OJ 



J? - k 



iri5(u>) 



TriS{u° - k) 



(B6) 



(B7) 
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together with Eq. 1|11[) . Comparing the left and right- 
hand sides of Eq. I|B6[) we see that the imaginary part 
should vanish, so we get 



1 + C0S<I>°X2l = 



(B8) 



which proves Eq. (|B1|I . In a similar way, starting from 
the equation for z a , 



Za=UJl + 2 Jo d \z a -k) + 



(l-cosfcx 2 i) (B9) 



we get 



1 



COS£l>°X 2 i 







(BIO) 



which proves Eq. I|B3|) . 

The uj° satisfy the integral equations 



cji + 2 / dk\ 2 vl- 



V 



■" - u! + 2 / dkX 2 v 2 k - 



u)°-k 
V 



Q°-k 



1 + cos 
1 — COS 



<2n + lWfc 



2mrk 



(Bll) 



Using graphical methods it can be shown the first equa- 
tion has a unique solution for each integer n, provided 
that 



lo\ — 2 lim 



dkX v] 



.V 



1 + cos 



(2n + l)7rfe 



> 
(B12) 



In the limit ui° — > the cosine term gives a vanishing 
integration. Thus Eq. (|B12|) is satisfied if Eq. is sat- 
isfied. A similar argument applies to the second equation 
in (Effi) . 

Eqs. IjBljl and l|B3|l explain the oscillatory behavior of 
7j seen in Fig. |S] 

To explain the oscillations of tUj we note that the terms 



inside brackets in Eq. (|B11|) are even in k around u>°, re- 
gardless of n. On the other hand, the principal parts 
are odd. Hence the product is odd and the integration 
around ui° vanishes. Thus the largest contributions to the 
integrals come from the tails of the principal parts. The 
"1" terms inside the brackets give a much larger contri- 
bution than the "cos" terms, because the latter oscillate 
with k. Neglecting the "cos" terms we get 



(B13) 



where u>\ is the one-atom shifted energy (see Eq. (JHJl). 
This shows that the ujj have approximately the same val- 
ues when their respective 7j vanish. From Eq. I|B13() we 



conclude that the period of the oscillations of 7j and ujj 
is approximately 2iz/u)\. 

Adding Eqs. (|B5J| and l)B9|l we see that for weak cou- 
pling the poles of the one and two-atom Green functions 
obey the relations 



Zl « (B14) 

So both cjj and 7j oscillate around the one-atom Cj\ and 
7i, respectively. 

Finally, we show that the "force" T s between the atoms 
is a maximum when the decay rate is zero, as seen in Fig. 
IT31 When 7 S = 0, we have 



dT° s 


d 2 £o° 


dx 2 \ 


dx 2 i 






2 


dk\ 2 vl- 


Jo 





, 9 (2n+l)7rfc . 
-fc 2 cos^ -r- 2 — (B15) 



As argued above Eq. (|B13p the integral of the cosine is 
small. Hence we have 



(B16) 



dT° 

dX21 

A similar argument may be applied to T a . 



APPENDIX C: SUB-RADIANCE IN d > 1 
DIMENSIONS 



In one dimension, the vanishing of the collective decay 
rate occurs for distances given by the conditions 



1 + (Tj cosd>°a;2i = 



(CI) 



Assuming the potential v k is rotationally invariant, in 
d > 1 dimensions, analogous conditions would be 



Q(d)dd[l + (Tj cos(w°a;2i cos 6)] = (C2) 



where 9 is the angle of the wave vector k, with respect 
to the line joining the two atoms. The function O is 2 for 
d = 1 and sin 9 for d = 3. 

We see that Eq. l|C2ll can only be satisfied for the an- 
tisymmetric state with Uj = —1 and for short distances 
X21 -C ujj ~ u>i 1 . This agrees with the results of Stephen 
anticipated by Dicke |l|. 
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